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Abstract 

We report on the detection of SiO masers in Asymptotic Giant Branch variables toward bulge/disk glob- 
ular clusters. In five out of six cases, the radial velocities are compatible with the optically measured radial 
velocities of globular clusters in the assessed uncertainty. Two sources, toward Terzan 5 and Terzan 12, lie 
very close to the cluster centers. The objects toward Pal 6 and Terzan 12 have luminosities appropriate to 
the AGB tip in globular clusters, while those toward NGC 6171, Pal 10, and Terzan 5 are brighter than 
expected. It is suggested that the latter three may have evolved from merged binaries, offering a test for 
binary-evolution scenarios in globular clusters, if the membership is approved. 

Key words: circumstcllar matter — globular clusters: general — radio lines: stars — stars: AGB and 
post-AGB 



1. Introduction 

Globular clusters (GCs, hereafter) are well-studied stel- 
lar systems containing the oldest stars in the Galaxy, and 
much of our knowledge of stellar dynamics and stellar evo- 
lution rests on them. With an age of about 10 Gyr, 
low-mass 1 Mq) stars are now on the AGB phase 
in GCs. While younger, more massive, AGB stars of- 
ten turn into infrared stars with thick dust shells due 
to heavy mass loss, the low-mass stars found in GCs do 
not appear to evolve in this manner. Although several 
red giants in GCs show small mid-infrared (MIR) ex- 
cesses that indicate weak mass loss (Origlia et al. 2002; 
Ramdani, Jorissen 2001), they bear little resemblance to 
the dust-enshrouded stars found in young Magellanic clus- 
ters (Tanabe et al. 1997). Maser emission is another char- 
acteristic often seen in such objects (see Habing 1996). 
However, no maser sources have ever been found in GCs 
(Knapp, Kerr 1973; Bowers et al. 1979; Cohen, Malkan 
1979; Dickey, Malkan 1980). A possible exception is the 
OH maser source V720 Oph in NGC 6171 (Frail, Beasley 
1994). Feast and Whitelock (1994) suggested two possi- 
bilities for this object: to be a foreground star or to be an 
evolutionary product of a merged binary. Whereas previ- 
ous surveys searched for OH and/or H2O masers, no SiO 
maser detection has yet been reported (for negative re- 



sults with the VLA, see the Web page In this paper, 
we report on a result of SiO maser detections in AGB 
variables, possibly associated with GCs. 

2. Observations and Results 

Targets were selected from our program of near in- 
frared (NIR) observations of more than 140 GCs taken 
from the known galactic total of 150. NIR monitoring 
observations started in 2002 March using an array cam- 
era (SIRIUS) attached to the 1.4 m IRSF telescope at the 
South African Astronomical Observatory [see Nagashima 
et al. (1999) and Nagayama et al. (2003) for instrumental 
details]. During the course of this work, we found sev- 
eral new red variables with (J — K)o > 3.0 mag. These 
stars are redder than the samples with small infrared ex- 
cesses, studied by Origlia et al. (2002) or Ramdani and 
Jorissen (2001), at least in (J — A^)o, so they are more 
promising candidates for SiO maser emission. 

For thirty two stars, the SiO maser lines (J = 1-0,?; = 1 
and 2) at 43.122 and 42.821 GHz (Lovas 1992), respec- 
tively, were searched for with the Nobeyama 45 m radio 
telescope during 2004 March 28-30 and on May 8. The 
half-power beam width was about 40" at 43 GHz. We used 
a cooled SIS receiver (Tgys ^ 180 - 250 K) and acousto- 
optical spectrometers with high (40 kHz; AOS-H) and low 

^ (http:/ /www. computing. edu.au/~bvk/astronomy/HET608/ 
project/). 
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Fig. 1. Spectra of the SiO J = 1-0, f = 2 (top) and v = l line (bottom). The observation dates are listed in the format YYMMDD.D 
for all objects. 



No. ] 

(250 kHz; AOS-W) resolutions having 2048 channels each. 
The spectrometer arrays covered velocity ranges of ±390 
km s-i and ±800 km s'^ in AOS-H and AOS-W, respec- 
tively. The conversion factor of the antenna temperature 
to the flux density was ^ 2.9 Jy . 

We detected SiO maser emissions from six red variables 
toward six GCs (among 32 surveyed stars) , mostly in the 
direction of the galactic bulge. For the previously known 
variable stars in NGC 6171 and Pal 10, the variable IDs 
(VI and V2, respectively) were taken from the catalog of 
variable stars in GCs by Clement et al. (2001). NGC 6171 
VI is V720 Oph, mentioned in section 1. The other four 
variables were newly found from our NIR observations, 
and named according to Clement's numbering system 
(namely Lillcr 1 VI, Pal 6 VI, Terzan 5 V5, and Terzan 12 
VI). Lines of both v = 1 and v = 2 were detected, except 
in the case of Terzan 12 VI, for which only v = 1 was 
seen. Their spectra from AOS-H are shown in figure 1. 
The detection criteria used by Izumiura et al. (1999) were 
satisfied by these lines, except for marginal detections of 
the w = 2 fines in Lillcr 1 VI and NGC 6171 VI. The 
integration time for each object was about 30 minutes, 
similar to those made in the past SiO surveys (Izumiura 
et al. 1999; Deguchi et al. 2000ab). The detections were 
confirmed in AOS-W spectra in all cases. The observa- 
tional properties of the detected objects are summarized 
in table 1, which contains the source positions, angular 
separations from the cluster center (rsio), SiO transitions, 
radial velocities (T^lsr), peak antenna temperatures (Ta), 
integrated intensities, and rms noise levels for both transi- 
tions from the AOS-H spectra. Table 2 lists the properties 
of the GCs. They were mostly taken from a compilation 
by Harris (1996). In the table, / and b arc the galactic co- 
ordinates of the cluster center, is the cluster velocity 
with respect to the Local Standard of Rest (with the es- 
cape velocity Vh at half radius between parentheses),^ 
and rt are the half-mass radius and the tidal radius of the 
cluster respectively, and [Fe/H] is the metallicity. In addi- 
tion, we give DM (the distance modulus of each clusters) 
and E{B — V) (the reddening) used to calculate the ab- 
solute bolometric magnitudes in subsection 3.3, and N2M 
and -Pcontm I the number of NIR SiO candidates within 34' 
of the cluster and the contamination probability, which are 
discussed in subsection 3.2. 

3. Discussions 

3.1. Position and Velocity 

All six sources are located well within the projected 
tidal radii, rt of the GCs. Terzan 5 V5 and Terzan 12 VI 
lie within the half-mass radii, rh, i.e., very close to the 
cluster centers within 1'. 

Because the radial velocity of SiO maser emission is 
known to coincide with that of the central star within a 
few km s~^ (Jewell et al. 1991), it can be compared di- 
rectly with the velocity of the cluster. Except for Lillcr 1 

^ Taken from the table of escape velocities used for 
Gncdin et al. (2002), which is available at http://www- 
astronomy.mps.ohio-state.edu/ ~ognedin/gc/ . 
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VI, the SiO radial velocities fall near to the optically mea- 
sured velocities of GCs with velocity separations of 8 - 37 
km s~^. The velocity dispersion within a GC depends 
on its mass; about 10 km s~^ for the most massive clus- 
ters, and only a few km s~^ for small ones (see table 2 of 
Mandushev et al. 1991). Gnedin et al. (2002) calculated 
the escape velocities (as well as the velocity dispersions) 
at half radii and the cluster centers of 147 GCs, using 
isotropic King models with a constant mass-luminosity 
ratio (M/Lv = 3); these escape velocities at the half radii 
are listed between the parentheses in the 4th columns of 
table 2. For NGC 6171 VI and Terzan 5 V5, the ob- 
served velocity differences of the SiO masers to the GCs 
are well within the escape velocities of the globular clus- 
ters, though they are slightly above the one-dimensional 
velocity dispersions given in the Gnedin's table. For Pal 6 
VI and Pal 10 V2, the velocity differences are a factor of 
about 1.8 larger than the escape velocity, and for Terzan 
12 VI, a factor of 4.2 larger. However, note that the ac- 
curacy of the radial-velocity measurements is quite poor 
for these three clusters. They were obtained from the av- 
erage of the radial velocities of a few giants, which were 
selected from a sample of a handful of stars, excluding 
objects with a large velocity separation; historic exam- 
ples in velocity measurements of GCs suggest a typical 
uncertainty of about ±20 km s~^ (sec table 3 of Coelho 
et al. 2001). For Pal 6, the velocity was derived from only 
two (table 3 of Minniti et al. 1995) or four giants (Lee 
et al. 2004); their values differ by about 20 km s^^. For 
Pal 10 and Terzan 12, the radial velocities of only two 
stars were used (Cote 1999). 

With this caveat, the velocities of the SiO masers are 
considered not to be incompatible with those of the clus- 
ters, except in the case of Lillcr 1 VI, whose membership 
is rejected. Better optical velocity measurements of the 
surrounding giants would be necessary to establish the 
memberships of the SiO masers from the velocity coin- 
cidence with greater certainty. On the other hand, the 
radial velocity of NGC 6171 has been studied based on 
the spectra of a number of giants (Da Costa, Seitzer 1989; 
Suntzeff et al. 1993; Piatek et al. 1994). Harris (1996) 
gives a velocity of —20.4 km s~^, obtained by combining 
their results. This is in reasonable agreement with the SiO 
maser velocity (—28.4 km s~^), despite the large angular 
separation (~ 9') from the cluster center. There have also 
been proper-motion studies concerning this cluster, which 
support the membership (Frail, Beasley 1994). Also, for 
Terzan 5, the radial velocity, which was recently derived 
using 4 stars (Origlia, Rich 2004), coincides with the pre- 
vious value listed by Harris (1996), as well as with the SiO 
maser velocity within a tolerable range. Thus, NGC 6171 
VI and Terzan 5 V5 cannot be rejected as being members 
of the clusters from a kinematic point of view. 

3.2. Concentration of Detections to the Galactic Bulge 

Among the six GCs, five are located in the galactic 
bulge, the exception being Pal 10, which is in the galactic 
disk. Because of the high density of AGB stars in these ar- 
eas, possible contamination by the bulge AGB stars with 



SiO Maser Sources toward Globular Clusters 



4 



Matsunaga et al. 

Table 1. SiO mascr properties. 



Object 


RA(J2000) 


Dcc(J2000) 




Transition 




Peak Ta 


Intcg. int. 


rms 




h m s 


f It 


(') 







(km s-^) 


(K) 


(K km s-^) 


(K) 


LiUer 1 VI 


17 33 35.2 


-33 22 01.8 


2.59 


1-0, 


1 


-275.1 


0.143 


0.337 


0.039 










1-0, 


2 


-275.1 


0.162 


0.126 


0.049 


NGC 6171 VI 


16 32 24.6 


-13 12 01.3 


8.98 


1-0, 


1 


-28.1 


0.185 


0.224 


0.035 










1-0, 


2 


-28.8 


0.230 


0.201 


0.048 


Pal 6 VI 


17 43 49.5 


-26 15 27.9 


2.67 


1-0, 


1 


156.8 


0.253 


0.604 


0.049 










1-0, 


2 


157.2 


0.229 


0.739 


0.045 


Pal 10 V2 


19 17 51.5 


+ 18 34 12.6 


2.51 


1-0, 


1 


12.7 


0.193 


0.332 


0.033 










1-0, 


2 


14.7 


0.173 


0.653 


0.041 


Terzan 5 V5 


17 48 03.4 


-24 46 42.0 


0.34 


1-0, 


1 


-106.8 


0.232 


0.778 


0.043 










1-0, 


2 


-105.9 


0.246 


0.868 


0.051 


Tcrzan 12 VI 


18 12 14.2 


-22 43 58.8 


0.64 


1-0, 


1 


76.9 


0.208 


0.621 


0.045 










1-0, 


2 








0.065 



Table 2. Cluster properties. 



Cluster 


I 


b 




rii 




DM 


E{B-V) 


[Fe/H] 


^2m(^MSx) 


p<I'SiO 

contm 




n 


n 


(km s"^) 


(') 


(') 


(mag) 


(mag) 








LiUcr 1 


354.84 


-0.16 


61.0 (19.5) 


0.45 


12.57 


14.81 


3.06 


0.22 


80(187) 


0.207 


NGC 6171 


3.37 


23.01 


-20.4 (17.2) 


2.70 


17.44 


14.02 


0.33 


-1.04 


l(5)t 


0.034 


Pal 6 


2.09 


1.78 


193.7*(20.9) 


1.06 


8.36 


13.82 


1.46 


-1.09 


52(96) 


0.148 


Pal 10 


52.44 


2.72 


-13.3*(13.6) 


0.99 


3.08 


13.79 


1.66 


-0.10 


5(19) 


0.014 


Terzan 5 


3.84 


1.69 


-82.4 (32.2) 


0.83 


13.27 


14.97 


2.15 


0.00 


64(70) 


0.003 


Terzan 12 


8.36 


-2.10 


106.3* (7.8) 


0.84 


3.10 


13.33 


2.06 


-0.50 


34(64) 


0.006 



* The value is assessed to be uncertain more than 10 km s ^ (see subsection 3.1). 

t Evaluated from IRAS PSC and 2MASS with A' < 10 and iJ - ii' > 0.6 (MSX not available at this b). 



SiO masers in the line of sight must be examined. The 
surface number density of SiO maser sources is estimated 
to be approximately 10-20 sources per square degree in 
the galactic bulge; this number was obtained by averag- 
ing the numbers found in the inner bulge in past sur- 
veys covering an area of 25° x 6° (Deguchi et al. 2000ab), 
corrected (with multiplying by a factor of ten) for lim- 
ited color ranges [0 < log(_F25/^'i2) < 0.1] and flux den- 
sity limits (Fi2 > 2 Jy) of the above-noted observations. 
This gives 0.16 at most for the number of spurious SiO 
sources falling within 3' from the center of a cluster in the 
bulge. Furthermore, a probability of contamination by the 
disk/bulge objects was evaluated for individual directions 
of GCs by counting MSX sources with 12 /xm flux den- 
sity above 1 Jy having a NIR counterpart falling within 
5" with K < 10 and H-K> 0.8 in the 2MASS catalog; 
these are potential candidate SiO maser sources, and from 
past experience, half of these objects exhibit SiO masers. 
The results are shown in the last two columns of table 
2, listing the number of NIR counterparts satisfying the 
above conditions within 34' (1 square degree) of the GC 
(between parentheses, the number of MSX sources above 
1 Jy in the same area), and the probability of contam- 
ination, finding more than one such object within rgio 
toward the globular cluster with a SiO detection rate of 
50%. These numbers are quite consistent with the value 
obtained from the above estimate based on the known SiO 



maser sources in the bulge. 

In addition, the radial velocities of the bulge SiO maser 
sources spread over approximately ±200 km s^^. The 
chance of an accidental coincidence of a bulge SiO source 
both in position and velocity is thus considered to be quite 
small unless the motions of the bulge GCs are dominated 
by bulge-bar star kinematics, though this is somewhat du- 
bious; e.g., Cote (1999). While that in Liller 1 is probably 
in fact an accidental coincidence, it is highly unlikely that 
all of our detections were by chance. 

It is also interesting to note that the present detections 
of SiO masers are limited to the bulge/disk GCs. As men- 
tioned in section 1, a previous massive SiO maser search 
with the VLA failed to find SiO maser emission in any 
halo GCs. Furthermore, mass-losing AGB candidates, i.e., 
the very red stars found in our NIR observations and/or 
those identified in the MSX/IRAS surveys, are only found 
in bulge clusters. These facts indicate that the bulge GCs 
have AGB stellar populations that differ in their charac- 
teristics from those in the halo. A possible explanation 
is that the bulge GCs are relatively metal- rich (~ 1/10 
to 1 solar abundance; see table 2). The molecular abun- 
dances, especially for SiO, are enhanced in a metal- rich 
environment. 
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3.3. Bolometric Magnitude 

Because the peaks of the spectral energy distribution 
of mass-losing AGB stars, such as those discussed here, 
lie in the region of the NIR and MIR, their bolomet- 
ric magnitudes can be estimated by combining our NIR 
photometric data with MIR catalogs, such as IRAS and 
MSX (Egan et al. 2003). For the newly discovered vari- 
ables, we adopted mean magnitudes from our NIR light 
curves. They all show large amplitudes that are typical 
for mass-losing stars. For the other two variables that are 
too bright for us to monitor, we adopted their magnitudes 
from the 2MASS AU Sky Data Release (Cutri et al. 2003). 
Comparing the magnitude of NGC 6171 VI with that re- 
ported by Feast and Whitelock (1994), it is probably not 
far from its mean magnitude. The adopted magnitudes 
are listed in table 3, which also indicates the correspond- 
ing MSX source and/or IRAS point source, if any. The 
bolometric magnitudes, mboi were calculated by integrat- 
ing under trapezoidal lines drawn in Xogv-vFn diagrams. 
While we ignored the contributions from shorter wave- 
lengths, we applied Ray leigh- Jeans curves for the longer 
ones. The uncertainties arc expected to be about 0.2 
mag, mainly arising from the uncertainties of the mean- 
magnitude determinations in the NIR and the flux density 
errors in the MIR. In table 3 the results are also listed 
along with the derived absolute bolometric magnitudes 
obtained by correcting for the assumed cluster distance 
moduli and reddening values. 

The derived Mboi values should be lower than those 
of the AGB tips, which depend on the age and chem- 
ical composition. Aaronson and Mould (1985) showed 
that the younger AGB stars evolve upwards in luminos- 
ity. They consider that the tip luminosity for metal-rich 
GCs is about —4.5 mag at an age of about 10 Gyr. Only 
Pal 6 VI and Terzan 12 have a Mboi value consistent 
with this limit, i.e., lower than —4.5. On the other hand, 
the merger of binary systems or direct collision can pro- 
duce massive objects that would reach higher luminosity 
(Bailyn 1995). This possibility is also mentioned by Feast 
and Whitelock (1994). This is related to a problem con- 
cerning planetary nebulae (PNe) in GCs. So far, four PNe 
are known to exist in GCs (Ps-1 in M 15, Pease 1928; 
GJJC-1 in M 22, Gillctt et al. 1989; JaFu-1 in Pal 6 and 
JaFu-2 in NGC 6441, Jacoby et al. 1997). While this num- 
ber is less than expected based on the assumption that ev- 
ery dying star produces a PN, current models predict that 
the masses of white dwarfs in GCs are too small to produce 
PNe (see Jacoby et al. 1997; Alves et al. 2000, and refer- 
ences therein). These authors suggested that only spe- 
cial objects, whose masses are augmented in close-binary 
interaction, can become PNe. Because the variable stars 
with heavy mass-loss have a lifetime roughly equal to that 
of planetary nebulae, it seems natural to expect the same 
number of them, both of which are considered to be prod- 
ucts of merged binaries. Since blue stragglers are common 
in GCs (Piotto et al. 2004) and also with the recent discov- 
eries that they contain numerous low-mass X-ray binaries 
and millisecond pulsars (Lyne et al. 2000), it is strongly 



suggested that close-binary formations and merger events 
occur relatively often in them (Pooley et al. 2003). With 
these indications and our new detection of masers with 
high luminosities (NGC 6171 VI, Pal 10 V2, and Terzan 5 
V5), we support the idea of a merged binary by Feast and 
Whitelock (1994). It may be slightly premature, however, 
at this stage to discuss their statistics and formation rates 
before further complete surveys are made. 

4. Conclusion 

Among the six SiO maser sources that we detected in 
the present survey, five, with the exception of Filler 1 VI, 
stand as potential candidates for GC members. Among 
them, the variables toward NGC 6171 and Terzan 5 are 
highly likely GC members, though those toward Pal 6, Pal 
10, and Terzan 12, are somewhat debatable. The member- 
ship must be approved (or disproved) much more firmly 
in future observations, because some uncertainty remains 
as to the coincidence of the maser velocities with those of 
the clusters due to a lack of precision in the optical veloc- 
ity measurements of GCs. These maser sources are dust- 
enshrouded AGB stars that were not expected in GCs in 
the past. Some of them have luminosities that are too 
high to fit the standard evolutionary tracks of low-mass 
stars, and may have evolved from binary mergers. On the 
other hand, SiO sources toward Pal 6 and Terzan 12 have 
luminosities that are consistent with the normal stellar 
evolution of old stars, though they might also be special 
objects at a certain stage of close-binary interaction. 
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